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Abstract —A strong motivation of charging depleted battery 
can be an enabler for network capacity increase. In this light we 
propose a spatial attraction cellular network (SAN) consisting of 
macro cells overlaid with small cell base stations that wirelessly 
charge user batteries. Such a network makes battery depleting 
users move toward the vicinity of small cell base stations. With 
a line adjustment of charging power, this user spatial attraction 
(SA) improves in spectral efficiency as well as load balancing. We 
jointly optimize both enhancements thanks to SA, and derive the 
corresponding optimal charging power in a closed form by using 
a stochastic geometric approach. 

Index Terms —Wireless power transfer, spatial attraction, spec¬ 
tral efficiency, load balancing, rate coverage, stochastic geometry. 


I. Introduction 

With the mobile data traffic explosion, the small cell net¬ 
work has come to the fore to offload over-utilized macro cell 
traffic. However, a small cell base station (BS) frequently 
serves a limited number of users due to its small coverage 
area 0- 0- In addition, the mobility of users even make the 
time spent in small cell area relatively short. This results in 
under-utilization of small cell resources |2;|, motivating us to 
balance traffic between macro and small cells. 

As a remedy for this, we suggest to enable small cell BSs 
to charge mobile devices wirelessly, given strong desire of 
mobile users to charge their devices |3j. Seeking for charging 
depleted batteries even makes users move to a nearby place 
providing power outlets 0 - Therefore, a wireless power 
transfer (WPT) in small cell BSs can act as an incentive to 
yield the spatial attraction (SA) of users in macro cells toward 
the vicinity of small cell BSs. Throughout this paper, we call 
this network a spatial attraction cellular network (SAN) and 
its attraction amount is determined by the charging power. 

The physical displacement from macro to small cells can 
be utilized for the improvement in macro-to-small cell load 
balancing. From the small cell users perspective who were 
already in the small cell coverage (see Fig. 1), this attraction 
also leads to the improvement in spectral efficiency. This is 
because the attraction toward their associated BSs not only 



Fig. 1. Illustration of an SAN for different charging powers (or charging 
range r c ) with the maximum SA distance r a ■ A small cell BS B 2 spatially 
attracts battery depleting users Ua to its charging rim (yellow thick circle). 

increases their received powers of the desired signals but also 
coincidentally decreases interference. 

In an SAN, these load balanging and spectral efficiency im¬ 
provements vary with charging power, but their corresponding 
behaviors are non-trivial. From the load balanging perspective, 
increasing charging power for instance may incur too much 
associations at small cell BSs. To mitigate such an issue, it 
requires an optimal charging power that balances the associ¬ 
ations between macro and small cell BSs. From the spectral 
efficiency perspective, on the contrary, high charging power 
increases the SA users’ distances from the associated small cell 
BSs, at which the charged powers are still satisfactory in that 
distance (compare (a) and (b) in Fig. 1). As a result, increasing 
charging power decreases spectral efficiency improvement. In 
this paper we thus provide the SAN rate maximizing charging 
power that jointly optimize the load balanging and spectral 
efficiency improvements. 

Feasibility of SAN - SAN implementation hinges on 
whether SA is viable in practice. Although mobile device 
proliferation is relentless, improving battery capacity has still 
been a major technical bottleneck |6). It makes users more 
prone to be in danger of battery depletion, leading to a strong 
battery charging motivation in reality |3j|. Such a motivation 
results in SA toward charging places, for example, when 







determining where to go for coffee 0> 0 or even when 
charging is the sole reason for movements |5j|. 

Another key enabler is to implement WPT capable small 
cells. By means of the advanced signal processing technologies 
such as an adaptive beamforming, microwave WPT has a 
potential to provide battery charging for most of small cell 
coverages while guaranteeing user safety |[8) . When utilizing 
a beamforming antenna having an aperture of 3 m in radius, 
for instance, it theoretically yields charging at the distance of 
10 meters with the charging amount of 0.5 Watt when the 
transmit power is 50 Watt ©■ 

Related Works - WPT aided cellular networks have re¬ 
cently attracted attention, accompanied by rapid advancement 
in WPT technologies Its advantage in network 

capacity perspective has been explored in (TTJ, (Tg. The 
network capacity enhancement comes from allowing more 
communicating users in that WPT mitigates battery depletion. 
Along with such improvement, we also capture the network 
capacity increment via increasing spectral efficiency owing to 
SA. The SA also provides improvement via macro-to-small 
cell load balancing without degrading spectral efficiency. It 
is unlike preceding schemes such as cell range expansion 
(CRE) [131, [J4] that balance the load in return for the 
spectral efficiency decrement. In detail, offloaded users in CRE 
are associated with small BSs who give lower signal power 
than the originally associated macro BSs. Than the originally 
associated BSs turn into dominant interferers, whose powers 
are even higher than the recived signal powers, leading to the 
decrement in spectral efficiency. 

Contributions - This study examines the ramifications of 
WPT in downlink small cells in terms of maximizing rate- 
coverage, the proability of a typical user achieving data rate 
greater than a threshold. The main contributions are listed 
below. 

• We propose an SAN that spatially attracts users to 
small cell BSs by the aid of battery charging WPT and 
thereby elevates the rate coverage of the network. We 
use a stochastic geometric approach to derive its rate 
coverage and demonstrate that significant improvements 
are feasible only with short SA distance, e.g. 125 % rate 
coverage when maximum SA distance is only 1 m (see 
Propositions 1 and 2 as well as Fig. 3 in Section III). 

• An optimal charging power is derived in a closed form 
in an asymptotic case where user density is much higher 
than that of BSs such as an urban area, providing a design 
guideline of SAN (see Corollary 1 in Section III). 

• To analyze SAN rate coverate, we consider data usage 
pattern dependent on user’s residual energy. The analytic 
result shows that the more users’ downlink (DL) usages 
are affected by their residual energies, the more rate 
coverage SAN can achieve (see Fig. 4). 


TABLE I 

Summary of Notations 


Notation 

Meaning 

Ut (or U h ) 

Low (or high) battery level user 

Bk 

fc-th tier base station for k E {1,2} 

$k (or $„) 

Bk (or user) locations 

(or \ u ) 

Bk (or user) density 

Pk 

B^ s information transmission power 

Pt C (or P rc ) 

Charging transmission (or reception) power 

Pu 

User energy consumption during a unit time slot 

r c 

Maximum charging range to charge P u 

r a (or f a ) 

User spatial attraction (or maximum) distance 

K (°r K) 

Rate coverage of a low (or high) battery user when 
associated with fc-th tier base station 

n 

Average rate coverage of user 


II. SAN System Model 
A. Information Transmission 

Consider a two-tier downlink cellular network comprising 
macro BSs B\ and small BSs B->- Let the subscript k G {1, 2} 
hereafter denote the corresponding tier. The fc-th tier BS 
coordinates follow a homogeneous Poisson point process 
(HPPP) with density \ k where <!>! and $2 are independent 
of each other. User locations <1> U also follow a HPPP with 
density A u , independent of $1 and <1>2. For simplicity without 
loss of generality, we consider a unity time slot between 
the realizations of the spatial locations. All notations are 
summarized in Table I. 

A BS in the fc-th tier If,- transmits a downlink information 
signal with power If. The transmitted signal experiences path 
loss attenuation with the exponent a > 2 and Rayleigh fading 
with unity mean. Both B 1 and B > utilize the same frequency, 
of the bandwidth W, yielding in inter-tier interference. Con¬ 
sider an interference-limited environment. By using Slyvnak’s 
theorem 115], the signal-to-interference ratio (SIR) at a typical 
user associated with a If; and located at the origin is given as 

.. ,,^1- 

SIRfc := 


P k h 


Bi 


r;=iE, 


.(O 


p./iW 
e-C n I 


B 


(0 r“ 


( 1 ) 


where f?[ 0) denotes the coordinates of the BS associated with 
the typical user, B^ for i G {1,2,---} the i-th nearest 
interfering BS coordinates from the typical user, and and 
li f ' ! the corresponding channel fading gains independently fol¬ 
lowing exponential distribution with unity mean. The notation 
| • | indicates the Euclidian norm. Multiple users in a BS are 
served according to a uniformly random scheduler CD- 


B. User Energy Consumption 

Users download information with probability u. Each in¬ 
formation reception consumes energy amount P u per unity 
time slot. For a user with residual energy amount L , it 
allows at most i times consecutive information receptions 











when iP u < L < (i + 1 )P U for an integer i > 0. We 
discretize L into T,’s according to the maximum number of 
information receptions i in a way that i = {-p-\ where |_-J is a 
floor function. For more clarity without loss of generality, we 
consider battery capacity is 3 P u , i.e. L £ {T 0 , L k , T 2 , T 3 }. 

Battery depleting users are assumed to attempt less frequent 
information receptions. This leads to specify the information 
download probability u as 


! 0 if L £ L 0 (Empty) 

up if L £ Li (Low) (2) 

u h if L£ {T 2j T 3 } (High) 


where Uh > up > 0. For notational brevity, let Up denotes 
battery depleting users whose residual battery energy L £ L\. 
Similarly, Uh indicates the users having L £ {T 2 ,T 3 }. 


C. WPT Charging 

Small cell BS T? 2 ’s are equipped with additional beamform¬ 
ing antennas dedicated for charging battery via microwave 
WPT. For charging battery, each user has an antenna which 
is associated with one of the predefined antenna patterns so 
that the main lobe point heads for her. Each B 2 transmits 
a charging signal with power P tc through a power transfer 
dedicated frequency whose bandwidth is assumed to be unity 
without loss of generality. The charging channel is separated, 
yielding no mutual interference. The transmitted charging 
signal experiences path loss attenuation with the exponent 
(3 > 2. Users associate with the nearest /i 2 . The received 
charging power is 

P IC :=G ro P tc max{|4 0) |,l}” /3 , (3) 

where G m denotes the main lobe gain. By means of highly 
sophisticated beamforming technique, the charging power re¬ 
ceived from non-associated /i 2 ’s is assumed to be negligibly 
small. For more brevity, we henceforth consider P lc is normal¬ 
ized by P u and always bigger than P u . Since residual battery 
level is elevated when P IC > P u , the maximum charging range 
r c becomes (. Pi c /P u ) 1 ^ which makes P IC be equivalent to P u . 

D. Spatial Attraction 

All users tend to keep their positions, but the battery deplet¬ 
ing users are willing to move for charging a unit energy amout 
P u with the SA distance r a (< f a ). Specifically, assuming f? 2 ’s 
broadcast their locations, a user is spatially attracted toward 
its nearest B> if the following two SA conditions hold: (i) 
L £ {To, L i} and (ii) r a < r a ■ The former indicates users are 
battery depleting (or Up), and the latter represents SA distance 
should be no greater than the users’ maximum feasible SA 
distance. The required amount of energy to receive the B 2 
location information is assumed to be negligible; so that even 
To users can successfully receive it. Such SA affects user 


c <(3) 



Fig. 2. Four-level residual battery energy Markov chain where the state Li 
corresponds to the maximum i number of information receptions. 


locations, and thus results in user mobility model specified 
by the following two phases. For each unity time slot, they 
occur in order. 

1. (Uniform Distribution): Users are uniformly distributed. 

2. (Spatial Attraction): SA users move toward their nearest 
/j 2 through the shortest paths until reaching Bf s charg¬ 
ing rim providing P u (see yellow thick circle in Fig. 1); 
non-SA users keep their locations. 


III. SAN Rate Coverage Maximization 


In this section we focus on WPT charging power P tc , as well 
as its corresponding r c , that maximizes average rate coverage 
1Z in an SAN, the probability that a typical user’s average 
downlink rate exceeds a target threshold 0. We rigorously 
specify the definition as 


TZ :— 


j-P( log[l + SIRfc] > 



(4) 


where N k represents the number of users in a If - associated 
with a typical user. 


A. Association Probability 

In SAN, a user is associated with a BS based on its residual 
energy and the received power strength. 


Lemma 1. The probability A^. (or A f k ) that Uh (or Up) 
associates with B k is given as 

-l 


ai = 


Afc 1 ^2 A; ( 


*=t 


P 


( 6 ) 


A{ = {-if [ A\e ] + {k - 1) (7) 


Proof. See Appendix-A. ■ 

Such a BS association changes the association distance as 
well as the number of associated users. They are derived by 
utilizing Lemma 1, to produce the rate coverage. 
















B. Rate Coverage and Its Optimal Charging Power 

To state our main result, the rate coverage, let q denote 
the steady-state probability vector of each residual battery 
energy level, q = [go, gi, < 72 , 93 ]> whose states and transition 
probabilities are shown as a Markov chain in Fig. 2. SAN rate 
coverage TZ is represented in the following proposition. 


Proposition 1. Rate coverage in an SAN is given as 


^ = E 



( 8 ) 


The residual battery energy dependent information reception 
probabilities Ph = (g 2 + g 3 )«/i and Pj J = q\ut> for gi, g 2 , and 
g 3 , which are obtained by solving qT = q where 


1-ELiC/W q(1) ce( 2) ce( 3) 

ui l-t^-X)i=i c *W c i(X) c ^( 2 ) 

0 ue l—ut—ci{l)ce{l) 

0 0 u h 1 -Uh 


Ce (2)=e-” X23 ~^ r * 2 

Cf (3)=l-e-’ rA23_? 

c h =l-e~^ 2 . 


_ 7rA 2 (r c +i ! a) 2 
c , 

g— ttA 2 2 ? r c 2 

rc2 , and 


The distribution of Bk associated user number is P(-/V/-=n)= 
Y2 P(Nj!=n— ?n)P(./V|=m) where 


m=0 


V / n!r(3.5) \ Afc 


h, \ n 
k 


K=*i)=- 


n 

r3.5 


3.5+ 


A u Ph^' “ (,l+4 - 5) 


\_3.5 3 - 5 r{nH.5) f X u PlAI 


n!T(3.5) 


A k 


3.5+ 


A k 

XuPlA^ - {n+i - 5) 


Given residual battery energy and associations, rate coverages 
are in 0 at the bottom of the page. 


Proof. See Appendix-B. 


It is worth noting that SAN rate coverage does not mono- 
tonically increase with the charging power P tc as shown in 
Fig. 3(b). The reason is explained by the following two 
perspectives. 

1) The trade-off between load balanging and spectral effi¬ 
ciency gains: Increasing P tc (or r c ) attracts more users 
toward small cell areas for load balanging while reducing 
the spectral efficiency improvement that can be achieved 
by shortening the distance to BS. 

2) .S’,4 motivation: Too small power P tc motivates few SA 
that cannot sufficiently operate an SAN so as to pro¬ 
vide a rate coverage improvement. Charging too much 
power, in contrast, retains mobile users’ residual energy 
enough and may lose their SA motivations, decreasing 
the improvements in both load balanging and spectral 
efficiency. 

This makes us turn our attention to derive the optimal 
charging power P t * maximizing rate coverage. This deriva¬ 
tion is not straightforward due to the user safety condition. 
The received charging power density of a user should not 
exceed the maximum safe power density ?/ ED- delimiting 
the feasible range of P t *. Furthermore, the maximum feasible 
amount of P t *, as well as its corresponding distance r c , is 
also restricted by P 2 Voronoi cell inscribed circle radius v\ 
otherwise, SA does not always increase spectral efficiency 
and load balancing, which may degrade rate coverage. We 
collectively consider such trade-off and constrains on P tc , and 
yield its optimal value. 


Proposition 2. Assuming the load at each BS is equal to its 
mean to simplify SAN rate coverage, the optimal charging 
power P t * is given as 


P t *=min< argmax^^ 

l P,c ' ' 


p H Ain't+p L Alnl 


^ i + 1 - 28 ^(p^+p^r 


p 7 ? 

r \C 5 -*tc 


(9) 


ni = 


n[ = - 


1 + P(g,M) 

-7tA 2 f-Ptc? +r a ) L 2 (0) 


-7tA 2 P tc P +r c 


nl = — 


U(9) A{ 

-ir\ 2 (pt c ? + r a ) L 2 (9) 


A% - Aie 


Ai 

2 

— 7T A 2 Ptc 
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o — 7rAir 2 [l+/)( 0)M )+Ap( 0 r )] 


dr 


1 I P (0,M) 
1_h A h 


— 7T A 2 Ptc ^ 


1 + 


P(0,M) 


— 7rA2 Ptc 


|i+ 


-r 

Ptc? 


U{0pAi 


A 2 A 2 + A 2 p(e,M) 


+ 


M 


6 — I" 00 f a \ 

P(e, x ) = (ew- 1)“ / # 2 (l + H 2 J 

■<( 7 ) ('»-!)"" V 


du,L k (0)='&$F&,M=(Pt c } +ra) (£)“ ,A=^ (^) ; 
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(a) For A u = 2xl0 4 users/km 2 , A 2 = 3xl0 2 BSs/km 2 (b) For A„ = lxl0 7 users/km 2 , A 2 =3xl0 3 BSs/km 2 

Fig. 3. Rate coverage in an SAN (Ai = 10 BSs/km 2 , Pi=43 dBm, P 2 = 23 dBm, W = 10 MFIz, 9=1 Mbps, a= 4, /3= 5, 97 = 10 W/m 2 , G m = 20 dBi, 
Gs = - 6.5 dBi (20)). 


where 


[!-e A27r -A 2 7rEi(-A 2 7r)], and 




/ 

r / P \ £i 

2 \ 

16A 2 +4 

Hi ). 

A,) 


( 10 ) 


where t rn denotes a beam width and Ei(x) an exponential 
integral function. 


Proof. Appendix-C ■ 

Although an optimal charging power P t * derivation is not 
possible in a closed form, it can be discovered through a 
linear search utilizing the above proposition, which is easily 
accessible compared to running an exhaustive simulation. 
Furthermore, to get a closed form optimal value, we consider 
an asymptotic case where the user density is much higher than 
BSs’. 


Corollary 1. For \ u ^>Ak, P£ is given as below. 
• When Ai -C A 2 , 


D* 

tn 


1 
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1 +P(0,T) 

A 

r a 

1 -| 

2 

1 A 

Gm 


07T 

Xl {^) a+X2 

12 
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2 r “ 
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• When Ai>A 2 , 


1 



^ 277 ^( 0 ^) 



+A 2 



The result provides a charging power control guideline in an 
SAN as follows. For low small cell density A 2 and/or the SA 
distance r a , the SAN is likely to suffer from macro cell traffic 
congestion, requiring to attract more users to small cells for 
load balanging. In such a scenario, P* c should be increased, 
and vice versa for the opposite situation. For a large threshold 
9, spectral efficiency dominates the rate coverage, and thus 
P t * should be decreased to shorten the association distances, 
increasing spectral efficiency. 

Fig. 3 illustrates the SAN rate coverage that is superior to 
the case with no SA, and the improvement increases along 
with SA motivation (or maximum SA distance f a ). Moreover, 
it captures the proposed network outperforms a network with 
CRE that is beneficial to load balancing yet is harmful to 
spectral efficiency G3- 11 4 ]. Specifically, when it comes 
to load balancing, CRE resorts to decreasing the received 
power from an associated BS. In an SAN, on the contrary, 
its load balancing coincidentally shortens user association 
distance, and therefore also increases the received power from 
an associated BS, leading to the further improvement in rate 
coverage than CRE. It is worth mentioning that if users are 
willing to move at least 20 cm for battery charging, an SAN 
become superior in rate coverage to the network with CRE. 
Though, SAN may achieve lower rate coverage than CRE if 
users are unlikely to be attracted, as shown in Fig. 3(a). 


Regarding a user safety requirement, the optimal charging 
power may exceed the maximum safe charging power (red 
dotted lines). However, such a situation can be detoured via 
high user density or deploying more small cells as in Fig. 3(b). 

Fig. 4 depicts the maximized SAN rate coverage mono- 
tonically increases with small cell BS density. This relentless 
improvement is expected to be different from the case with 


Proof. Appendix-D 


























— DL usage at high battery : usage at low battery = 1:5 



Fig. 4. Maximized rate coverage and optimal charging power with respect to 
the A 2 and DL usage ratio ^ (Ai = 10 BSs/km 2 , A u = lx 10 4 users/km 2 , 
Pi =43 dBm, P 2 =23 dBm?V=10 MHz, 6=1 Mbps, a = 4, (3= 5, r a = 
2 m, G m = 10 dBi, G S = -7A dBi |20]). 


Further extension to this work therefore should involve their 
joint optimization problem. Another interesting avenue for 
future work is a network economic analysis on the SAN. In 
a similar way to ED> the SAN rate coverage optimization 
will jointly incorporate BS density and spectrum amount for 
communications, also combined with such values for WPT, 
providing a network operator’s SAN deployment guideline. 

Appendix 

A. Proof of Lemma 1 

For Uh* the association probability A 1 /, is directly derived by 
using Lemma 1 in |jl9]. For Ug, the association probability A ( k 
should incorporate the SA impact. Let ?> denote the distance 
between a typical Ug and his nearest /i/ i; . Now that a Ug is 
attracted toward the nearest B 2 if r 2 <r s , A\ is given as 


CRE whose improvement is saturated for high small cell BS 
density where load balancing is not much needed (see the 
convergence of dotted black and solid green lines in Fig. 3(b)). 
SAN can increase rate coverage via its spectral efficiency 
improvement even viable in such an environment. Focusing 
on the exponentially decreasing optimal charging power along 
with small cell BS density, deploying small cells can be a 
viable solution for ever-growing rate coverage while abiding 
by a safety requirement. 

In addition, the proposed SAN can improve much more rate 
coverage when DL traffic usage ratio of low-to-high battery 
users is large, i.e. when users are more sensitive to their 
residual energy. This represents that SAN effects well along 
with a growth of booming wireless communication services 
that require large battery consumptions. 

IV. Conclusion 

In this paper we have analyzed an SAN that spatially attracts 
battery depleting users by means of providing WPT battery 
charging. The result sheds light on WPT charging power 
control in order to maximize rate coverage of the SAN. Our 
analytic result reveals that charging power increase as much 
as possible does not guarantees the maximum rate coverage, 
and thus its optimization is required. Such an optimal charging 
power is derived in a closed form, and its corresponding de¬ 
sign guideline is presented for different network deployment, 
wireless channel environment, and a user safety requirement. 

The weakness of this study is its one-sided focusing on 
a downlink scenario. SA has a potential to further improve 
uplink rate coverage by its shrinking association distances. In 
such a scenario, the uplink optimal charging power may not 
be in accordance with the downlink optimal charging power. 


A[=P(Ug —>Bi and r 2 >r s ) 


=P(P 1 r 1 a >P 2 r 2 a and r 2 >r s ) 

, re _m ' 2E “= lAi (^)“dr 


== 27tAi 


+27rAie~ 7rA2r 


, P 2 , 


re~* Xir2 dr (13) 


where (a) follows from the nearest neighbor distance distribu¬ 
tion 1151. Similarly, A 2 is given as 


A^=P(Ue-tB 2 or r 2 <r s ) 


=P (P 2 r 2 “<Pitt a , r 2 >r s )+P(r 2 <r s ) 
=2ttA 2 f re _7T ( E -= 2Ai ^ ' +X2 ) r2 dr+l-e-* 

J TVo 


(14) 


that finalizes the proof. ■ 

B. Proof of Proposition 1 

The desired result comes from the following four parts. 

1) Residual Energy Dependent Information Reception 
Probabilities Ph and Pry Let cg(i) denote the probability 
that Ug is capable of additional i times information receptions 
thanks to WPT battery charging. Its calculation is followed 

by the fact that Ug can receive i times information receptions 

_2 

thanks to battery charging if 2 £ r c <r 2 <r c +r a for i= 1 and 
if (i+2)~P r c <r 2 <(i+l)~ p r c for 1=2,3. Let Ch denote the 
battery charging probability of Uh ■ The event of c,/, is identical 
to the event that Uh is located with the distance closer than r c 
from the nearest B 2 . Applying void probabilities m to their 
spatial regions of cg(i) and Ch yield the calculations. 













2) BS Association Distance Distribution: Now that an SA 
user goes toward the nearest B 2 , the distance between them 
diminishes. For the SA affected association distance for Ug, 
its complementary cumulative distribution function (CCDF) is 

rot ^ \tt n \ P( r k>r and Ue->B k ) 
P(r k >r\Ug^B k )= - — t - (15) 

At 

where U^B k denotes an event that a user U associates with a 
B k . For the probability density function (PDF) (r), calcu¬ 
lating the numerator in CCDF requires a simple modification 
of ( fl4| in the proof of Lemma 1. For (r), the numerator 
calculation is divided into two cases: (i) if r c < r 2 < r s , an Ug 
is located at a distance of r c from P 2 ; (ii) if r 2 <r c , r c is 
restricted to the B 2 cell, making the calculation become the 
void probability JT5| of a ball with radius r 2 . Applying Lemma 
1 and differentiating CCDF derive fl{r) as below. 


2nXl re~ n ^= lXi (^) ‘ 
f*(r)={ re 

2ttXi c -TT\or, 2 rc ~TT\i r 2 


2tt\ 2 vr —n\or 2 


fvAr)={ 6 r 


-7r\ 2 r 


2n\2 

Ai 


re 


Ai, 


if r<r 
otherwise 


if r<r c 
if r c <r<r s 

otherwise 


(a) 4 . 


(16) 


(17) 


where S rc is a Dirac delta function yielding 1 for r=r c , 
otherwise 0. The PDF (r) for U k is provided by Lemma 4 


2 ^ 


in [141 as Ah 

J) Distribution of B k Associated User Number: The proba¬ 
bility mass function of the user number N k excluding a typical 
user associated with a B k is derived by using Lemma 1 and 
Corollary 1 in [14 j. 

4) Rate Coverages Conditioned on Residual Energy and an 


Association lZ k and lZ f k : Applying Theorem 1 in [ 141 with the 


preceding results and Lemma 1 enables to calculate the rate 
coverages. When it comes to averaging interferer locations, 
it is worth noting that small cell interferer distances at a B\ 
associated typical Ug range from r 8 to 00 if r\<M\ otherwise, 
the interferer distances range from 0 to 00. ■ 


C. Proof of Proposition 2 

For analytical tractability, we approximate the number of 
the fc-th tier associated users as its average value as below 

HD- 


E[N k }= l ^(P H A h k +P L A{). (18) 

The feasible range of charging power is restricted by a safety 
requirement. Since a user recieves the charging power with 


main lobe gain from his associated B 2 and received power 
from other P 2 ’s is assumed to be negligibly small, the average 
received power density P rc is represented as 


p — 


G m P t c 


7rtan 2 (l| i )max| | ,l| 
‘2G rn A 2 P\ c 


tan 2 (^) Wo 


r 2 (^ 
2G m \2Ptc 
tan 2 (lg L ) 


re 


— A 2 


dr+ 


00 re -Tr\ 2 r 


dr 


[l-e- A2 "-A 2 7rEi(-A 2 7r)] 


where (a) follows from surface area of power transmission 
in main lobe [20], t rn denotes a beam width, and Ei(i) an 
exponential integral function. 

Applying the safety requirement P rc f jf results in I f .. 
Another factor delimiting the feasible range of charging power 
is Voronoi cell area, i.e. the maximum charging distance r c is 
no larger than the inradius of a B 2 cell area v. To simplify 
our exposition, we approximate v as its expected value E[z/] 
calculated by using m- 

EM = (16A 2 +4[1 + (P 1 /P 2 ) 1/q ] 2 A 1 )- 1 / 2 . (19) 

Applying the maximum charging range requirement Pf.= 
completes the proof. ■ 


D. Proof of Corollary 1 

When A„>A fe , n^P L A{{ni-n 1 i){P H A 1 i+P L Ai)~ 1 . 
Applying Taylor expansion with a linear interpolation 
simplifies K as (4f3r c +6f 2 r c 2 +4f Q r c 3 +f Q ) - 

(l+P(e :T )) (2f Q r c +ff) for Ai<CA 2 . In a similar way, the 
approximation of TZ for Ai>A 2 becomes 7rA 2 (2f a r c +f 2 ) — 
7 t 2 A 2 2 P( ^ ) (2f a r e 3 +r 2 r e 2 ). Differentiating the approxi¬ 
mated 1Z with respect to r c yields the optimal charging power. 
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